Phase diagram of the su(S) quantum spin tube 



J. de Gier, M. T. Batchelor and M. Maslen 
Department of Mathematics, School of Mathematical Sciences, 
Australian National University, Canberra ACT 0200, Australia 
(February 1, 2008) 

We calculate the phase diagram of an integrable anisotropic 3-leg quantum spin tube connected to 
the su(8) algebra. We find several quantum phase transitions for antiferromagnetic rung couplings. 
Their locations are calculated exactly from the Bethe Ansatz solution and we discuss the nature of 
each of the different phases. 
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I. INTRODUCTION 

The properties of n-\eg quantum spin ladders have attracted considerable recent attention from both theorists and 
experimentalists jjj. The main reasons for this interest are that, like their one-dimensional counterparts, ladders 
are amenable to approximation schemes such as bosonization, and to accurate numerical investigation by the DMRG 
method. It is now well established that the n-leg Heisenberg ladders have the rather striking property of being massive 
for n even and massless for n odd. A number of ladder-like compounds have been found. For example, the 2-leg 
ladders SYCU2O3 Q and C^Cs-ff^A^^C^ || exhibit a gap, while the 3-leg ladder SY2CU3O5 || is gapless. 

Although the Heisenberg ladders are not solvable in the sense of the Heisenberg spin chain, a number of solvable 
ladder models are known [j^-jOl. One 2-leg ladder model with fixed rung, leg and diagonal interactions is equivalent 
to the XXZ chain at A = — |and is thus massive [||J|]. Another 2-leg model has been constructed which is massless 
in the absence of rung interactions with a transition to a massive phase at a critical (non-zero) value of the rung 
coupling 0. This is in contrast with the 2-leg Heisenberg spin ladder for which the transition to the massive phase 
is at zero rung coupling. In the integrable model the Heisenberg rung interactions appear as chemical potentials 
which break the underlying su(4) symmetry. The phase diagram of this model has been established by means of the 
Bethe Ansatz solution 0. The integrable model has been generalised to an arbitrary number of legs in which the 
underlying su(2 n ) symmetry is broken by the Heisenberg rung interactions, which again appear as chemical potentials 
in the Bethe Ansatz solution ||] . Integrable ladder models with Hubbard || and t — J Jl0| interactions have also been 
recently found. Others have been constructed with o(2 n ) and sp(2 n ) symmetry [jnj. Some ladder models have also 
been constructed with dimer or matrix product ground states |12j| . 

In this paper we investigate the thermodynamic properties of the solvable 3-leg qnantum spin tube || . 

II. THE MODEL 

We consider an integrable 3-leg spin ladder, or spin tube for periodic boundary conditions. The spins along each 
leg and each rung have an isotropic Heisenberg interaction, with the introduction of many-body terms to retain 
integrability. As usual, it is convenient to express the Hamiltonian in terms of the Pauli matrices, 



1\ v (0 -i\ z I 
10' = \ i ' a = U -1 



(2.1) 



The Hamiltonian of our model is 
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The operators (cr 1 )^, (<J y )f^ and (cr 8 )^ act as the corresponding Pauli matrices on the (i, Z)th factor in the Hilbert 
space, 
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We take periodic boundary conditions along and also across the ladder. Now we note that h*f is just the permutator 



on the space Vi (8> Vj and thus that ( |2.3[ ) is simply the integrable isotropic permutation Hamiltonian corresponding to 
the su(8) algebra |Q. Since [h*j, Zi™ ns + h\ L "^} = 0, H defined in ( p^ ) is also integrable. We set J\ = J2 = J and 
J3 = J', so that we can go from the isotropic tube, J' = J, to the ladder, J' = 0. This model thus generalizes that 
given in || to anisotropic rung couplings. 

It is convenient to change to the basis where the square and the z-component of the total spin of a given triangle, 
S = erW + er( 2 ' + er^ 3 ) are diagonal. It follows that the eight states on a given triangle fall into a spin-| quadruplet 

and two spin-i doublets. This basis, in terms of the (<7*)W eigenvalues, is, 
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It is emphasized that the Hamiltonian (2.3) has the same form on this new basis. The rung Hamiltonian (2.4) now 
becomes diagonal and is given by /i rung = diag{— 3 J, —3 J, —2 J' — J, —2 J' — J, 0, 0, 0, 0}. 

As claimed above, H can be diagonalized using the Bethe Ansatz. The Bethe Ansatz equations are well known 
|l3f , and given by 
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where j = 1, . . . , M r with r — 2, ... ,7 and M 8 = 0. The eigenenergies of J2i=i h% f+i are gi ven by 
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It is important to note that since the permutator (2.3) is invariant under any ordering of the states ( |2.6| ), this result 
may be obtained using any choice of reference state (or pseudo- vacuum) |f2) and any assignment of Bethe Ansatz 
pseudo particles. For each cho ice however, one has to re-interpret the numbers M r in terms of those corresponding 
to the ordering chosen in ( |2.6|) . The rung Hamiltonian does alter with the choice of ordering, but the change is just 
a rearrangement of its eigenvalues along the diagonal. We use this property to our advantage by doing calculations 
with that choice of ordering for which the Bethe Ansatz reference state is closest to the true ground state of the 
system. Below we list the six different possibilities corresponding to the choice of ordering of the two doublets and 
the quadruplet. 

1. {|0>, |1), |2), |3>, |4), |5), |6), |7», E = E^ + 2(J - J')M 2 + (2J' + J)M 4 . 
The Bethe Ansatz reference state is 

L L 
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and the assignment of pseudo particles is such that M k = Er=fc-^j wnere N r is the number of states \r) 
occuring in an excited state. In the cases below the reference state is understood to be the tensor product of 
the first state listed. 

2. {|2>, |3>, |0), |1), |4), |5), |6), |7}}, E = E lc & + 2(J' - J)M 2 + 3JM 4 . 

Now the pseudo particle assignment is given by Mi = N + Ni + N 3 + Er=4 N r > M 2 — N + Aq + Er=4 N rt 
M 3 = Aq + N r and M k = J2l =k N r for k > 4. 

3. {|2>, |3), |4), |5), |6), |7), |0), |1}}, E = E lc % + {2J> + J)M 2 - 3JM 6 , 

with M k = Er=o N r + EL/c+2 N r ^ 1 < * < 5, M 6 = N Q + N t and M 7 = Nl 

4. {|4), |5), |6), |7), |2), |3), |0), |1}}, E = E lc % - (2J' + J)M 4 + 2(J' - J)M 6 , 

with M k = Y? r =o N r + El= k +i N r for fc < 3 and M 4 = N + Nx + N 2 + N 3 , M 5 = N Q + N x + N 3 , M 6 = N Q + N 1 
and M 7 = Ny. 

5. {|4}, |5), |6), [7), |0), |1), |2), |3}}, E = E^ - 3JM 4 + 2(J - J')M 6 , 

with M k = Ei? =0 N r + Er=fe+4 for fc < 3 and M k = Er=fc-4 N r for * > 4. 

6. {|0}, |1), |4), |5), |6), |7>, |2), |3}}, E = E^ + 3JM 2 - (2 J' + J)M 6 , 

with Mi = JVi + M a , M fc = Er= 2 N r + Er=fc+ 2 N r 2 < k < 5, M 6 = N 2 + N 3 and M 7 = A^ 3 - 
Of course, all of these different choices for the reference state do not change the physics, but it turns out that each 



of them is convenient in some part of the phase diagram, as we will see in Section IV. Before calculating the complete 
phase diagram we first consider two interesting lines. 

III. SPECIAL LINES 
A. The tube: J' = J 

Due to the extra symmetry in this case, the two doublets become degenerate. Consider now the Bethe Ansatz 
where the reference state is the spin up state of either of the spin-i doublets. According to case 1. or 2. in Section 
O, the energy is given by 

Mi 

E = -J2-Jru- T+3./M 4 , (3.1) 

where M 4 denotes the number of quadruplet excitations. It is obvious that if J is large and positive, no quadruplet 
state can exist in the ground state. There are three Fermi seas and therefore three gapless excitations. All flavours 
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of the quadruplet excitations are massive and their gap may be calculated by creating a A^ 4 ) mode and a A( 3 ) hole in 
the Bethe Ansatz equation jjj (see Appendix The gap is 

A = 3J-~(log4 + 7r). (3.2) 

Thus we find a phase transiton from J > J+ = A (log 4 + tt) to < J < J+ where all quadruplet flavours become 
massless and there are seven Fermi seas and no massive excitations. Note that only at J = 0, all Fermi seas are 
completely filled. 

Considering the case J < 0, it is convenient to choose either case 4. or 5. of Section || The energy is given by 

Mi 

i=l l A j J + 4 

where now M4 denotes the number of excitations on the spin-| ground state. These excitations are the four flavours 
belonging to the two spin-i states and are all massive for J <C 0. We can use the same calculation as above to 
conclude that again there is a phase transition at J = J c = — where the two spin-^ doublet excitations become 
massless. 

B. The ladder: J' = 

As above, we start by considering this case for J > 0. It is most convenient to take the Bethe Ansatz as in case 1. 
of Section |j| The energy is thus given by 



M- 

E = -J2 — (TTT^ 7 + 2 Jil/ 2 + JM 4- ( 3 - 4 ) 



For J ^> it is obvious that M2 = for the ground state. The ground state is thus formed by the {|0), |1)} doublet 
and there is one gapless excitation and six massive ones. As before we can calculate the gap by introducing a A^ 2 -* 
mode and a A*- 1 ) hole, giving rise to 

A = 2J-2 J y^— -duJ = 2(J-log2). (3.5) 

This gap thus closes at J c \ = log 2 and for J < J C) j there are three massless excitations and four massive ones 
corresponding to the triplet states. These become massless at a second critical point, J C) 2, whose value can not be 
calculated analytically. This point continuously connects to the critical point J+ on the line J' = J. 

IV. PHASES 

As already noted, the different choices of the reference state do not change the physics, but each of them is convenient 
in some part of the phase diagram. To appreciate this, it is helpful to divide the phase diagram into six regions defined 
by the lines J = 0, J' = J and 2 J' + J = as in Fig. [[]. In each of these regions, the corresponding choice of reference 
state of Section || is the most convenient. 

1. J > 0, J - J' > and 2J' + J > 0. 



This part of the phase diagram contains the ladder for J > 0. As in Section III B we choose to describe the 
spectrum as in case 1. of the list in Section [ll| For J ^> J' it is obvious that for the ground state one must have 
M2 = 0. The ground state thus is formed by the doublet {|0), |1)}, which is effectively a spin-i XXX chain. 
The system is critical with one massless excitation and six massive ones. The gap can again be calculated by 
creating a A^ 2 ^ mode and a A^ hole (M2 — * 1, Mi — * M\ — 1) in the Bethe Ansatz equations, and is given by 

A = 2(J- J') - 2 log 2 - (4- 1 ) 

Thus there is a critical line on which this gap vanishes, given by 
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J' c = J -log 2. 



(4.2) 



Above this line some of the excitations on the {|0), |1)} doublet become massless. We now investigate which 
ones. From the form of the energy it follows that if 2 J' + J is large, M4 = and o nly the {|2), |3)} doublet 
excitations have become massless. Thus for large 2 J' + J and above the critical line (4.2), the ground state is 
given by the two spin-i doublets and the quadruplet excitations are massive. There is another phase transition 
where al so t he quadruplet excitations become massless, i.e. M4 7^ for the ground state, and its location on 
the line (4.2) is given by 

2J' C + J = Q, or (J, J') = (| log 2, -f log 2). (4.3) 

Its location changes as M2 increases, which can be calculated using perturbative calculations (see Appendix^). 
It ultimately reaches the value of 

(J, J') = (i(fog4 + n), i(log4 + tt)), (4.4) 

on the line J' — J where three Fermi seas are completely filled instead of just one. 

2. J > 0, J - J' < and 2 J' + J > 0. 

By the same reasoning as above the ground state is again formed by a doublet for J' J, but now {|2), |3)} 
since we use a different reference state. The system again is critical with one massless excitation and six massive 
ones. The gap can be calculated similarly and the critical line on which this gap vanishes is given by 

J^ = J + log2. (4.5) 

Below this line the {|0), |1)} doublet excitations have become massless. There is another phase transition where 
also the quadruplet excitations become massless, i.e. M4 ^ for the ground state, and its location on the line 
(4.5) is given by 

(J, J') = (0, log 2). (4.6) 

Its location changes as M2 increases reaching the previously calculated value, 

(J, J') = (£(log4 + tt), ^(log4 + tt)), (4.7) 

on the line J' = J where three Fermi seas are completely filled instead of just one. 

3. J < 0, J - J' < and 2 J' + J > 0. 
Here we find the phase boundary, 

2J£ + J = 21og2, (4.8) 

where the spin-| quadruplet excitations on the {|2), |3)} doublet ground state become massless. The point on 
this line where also the spin-i doublet {|0), |1)} becomes massless is given by 

J = 0. (4.9) 

As before, its location changes, reaching 

(J,J') = (~^), a 6 « 0.5076..., (4.10) 

3 D 

on the line 2 J' + J = 0. 
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4. J < 0, J - J' < and 2 J' + J < 0. 

In this region the ground state for large and negative 2 J' + J is the quadruplet. It's smallest excitations, those 
of the {|2), 1 3)} doublet, become massless at 

2J' + J= -i(7r + log4). (4.11) 

On this line the other doublet becomes massless at 

(J, J') = (-£(tt + log 4), -^(tt + log 4)), (4.12) 

on the line J' = J. This point moves to 

(J,</') = (-y,y), (4.13) 

on the line 2 J' + J = 0. 

5. J < 0, J - J' > and 2 J' + J < 0. 

In this region the quadruplet is still the ground state for large and negative J, but now the lowest excitations 
are the {|0), |1)} doublet. They become massless at 

J=-^(7r + log4), (4.14) 

independent of J'. The other doublet becoms massless at 

(J, J') - (-i(7r + log 4), -i(7T + log 4)), (4.15) 

on the line J' = J and this point extends to 

(J,J') = (0,-^), (4.16) 

on the line J = 0. 

6. J > 0, J - J' > and 2 J' + J < 0. 

Here, the ground state for large J is the {|0}, |1}} doublet, and the smallest excitations are the quadruplet which 
become massless at 

J = flog 2, (4.17) 

independent of J'. The other, doublet, excitations become massless at 

(J, J') = (f log 2, -| log 2), (4.18) 

on the line 2 J' + J = 0. On the other phase boundary, J = 0, this point is 

(J,j') = (0,-y). (4.19) 

Our findings, along with the phase boundaries, are summarized in Fig. |l|. 
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FIG. 1. Phase diagram. The dashed lines are J' = J and 2 J' + J = which together with the y-axis divide phase space into 
six regions. The dotted lines are sketches of phase boundaries. 

V. CONCLUSION 

I n th is paper we have calculated the phase diagram of the integrable anisotropic 3-leg quantum spin tube defined 
in ( |2.2[ ). The model differs from the usual Heisenberg spin ladder in that it includes additional four- and six-body 
terms in the leg Hamiltonian The phase diagram, as calculated from the Bethe Ansatz solution, is given in Fig. 

[I]. The straight line phase boundaries can be calculated exactly as well as the critical points on the lines where the 
rung Hamiltonian becomes partially degenerate. For the location of the dotted phase boundaries one has to resort 
to perturbative calculations which show that they end parallel to the degeneracy lines. In Fig. |l] we have sketched 
the full form of these phase boundaries. It is to be hoped that their precise location may be accurately determined 
numerically. 

Since the model has an odd number of legs, it is critical in the entire phase space and has no true gap. It is however 
the simplest example with more than one quantum phase transition for J > and serves as an example for what to 
expect for higher leg ladders. In general, if the rung Hamiltonian has n different eigenvalues, one may expect n—1 
different phase transitions along a generic line between the phase with the lowest and highest total spin. At each of 
these one expects to find quantum critical behaviour, which will however diminish in magnitude. 

The 3-leg quantum spin tube is also integrable in the presence of a magnetic field. The corresponding 3-leg 
Heisenberg model has been shown to exhibit magnetization plateaus by means of renormalization group calculations 
based on the bosonized effective continuum theory |]l4| . The exact calculation of the magnetic behaviour of the present 
model will thus be of considerable interest. 
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APPENDIX A: GAP CALCULATION 



As is well known, the Bethe Ansatz equations (2.8) may be written as integral equations in the thermodynamic 
limit (L — > oo), 



01(A) 



(A)+ / ' a 2 (A-A')p (1) (A')dA'- / o x (A - A> (2) (A')c!A' 



(Al) 



pW(A)+ / " a2 (A — A')p^(A')dA' 



(r-1) 



Oil 



(A- AV~ 1) clA'+ / 01 (A- A> (r+1) (A')dA', ( A2 ) 



where /?' r ) are the distributions of Bethe Ansatz roots and holes and 



ft [ 

a " (A) = 2^*TW (A3) 



In the sequel we will make use of the following convention for the Fourier transform, 



m = ^J_ /(fc)e" ifc "dfc, f(k) = J />)e ito du;, (A4) 

and of the results, 

a n (k)e- ikuJ dk = e-^ 1 , for n > 0, (A5) 

^^^=l^ m m-^°Hi)), ^ (0) (^)= r ^- ( Ag ) 

As an example we give here the calculation of the gap for the tube when J' = J > 0. The ground state for J ^ 
consists of three completely filled Fermi seas. The Bethe Ansatz equations ( |Al| ) and QA2|) may therefore be solved 
with Fourier transforms. A quadruplet excitation is created with a A^ 4 ) mode and a A^hole. Following Wang M, 
we will denote the changes in the distribution functions as a consequence of this particle-hole pair creation by 8pv) . 
The Fourier transformed Bethe Ansatz equations then become (with some abuse of notation) 

= 5p {1 \lu)(1 + e-W) - 6pW(u)e-$M 

5p^(u)(l + e-W) = 5^(u)e~^ u \ + 5p^\cu)e-^ - -Le^^M 

2nL 

S P ^(lu)(1 + e-W) - J_ e -^"-\"\ + Sp (2 \cu)e-^ + J-e-^-iM, 
2nL 2nL 



from which we obtain 



1 e -* W " - e h\->\-* W » 

6pV{u,) = — —— — — (A7) 



2nL 4 cosh lo cosh - 2 



Thus it follows that the smallest gap (A 1 - 3 -* — > 00, A^ 4 ) — > 0) is equal to 

/OO 
e-sl^l^W^dw = 3J- i(log4 + 7r). (A8) 
-OO 
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APPENDIX B: PERTURB ATIVE CALCULATIONS 



Suppose we are in region I in the phase where both doublet excitations are massless, but close to the phase boundary 
J' = J — log2. The {|2), |3)} excitations have just become massless and M 2 /L and M 3 /L have a small but finite 
value. The X^ Fermi sea will be filled up to some finite momentum A while those of A*- 1 -* and X^ are completely 
filled. By creating a X^ excitation, we change A and the distribution functions. The Bethe Ansatz equations for this 
situation are given by 

/ a 2 (A-A')p (1) (A')dA' = ai (A)+ / _ ai (A-A')p (2) (A')dA', (Bl) 



OO 

A /*+oo r+oo 



j o (2) (A)+ / a 2 (A- A')p (2) (A')dA' = / a x (A - A')p (1) (A')dA' + / <n(A - A> (3) (A')dA', (B2) 



— A J — oo 

+00 r +A 



P (3) (A) + / °° a 2 (A - A')p (3) (A')dA' = f ai (A - A> (2) (A')dA' + ^1 j (B3 ) 



00 



with the conditions 

r +A T./r_ r +oo 



J P (2 HX)dX = ^, /_V)(A)dA = ^. (B4) 



-A 

These equations can be solved perturbatively in A which we will do up to 0(A 3 ). With the help of 

,+A 



-A 



a n (X - A')p (2) (A')dA' = 2A a „(A)p( 2 '(0) + ±A 3 «(A)p( 2 ) ( ) + a „(A)p( 2 )"(0)), 

= a„(A)^ + iA 3 a^(A)p^(0), (B5) 



we find 



2 



p (2) (0 ) = 3°|2(i + I) + 0(A). (B7) 



A = A(i-i) + of4)- fBs < 



Since we want to keep M 2 constant it follows from (B7) that 
Finally, the energy is given by 

/+00 
j o (1) (w)e-5l"ldw + 2(J - J')M 2 + 2 J' + J 
-OO 

= -2Llog2(l + ^) + LA 3 C(3)1 ° g2 (l + I) + 2(J - J')A/ 2 + 2 J' + J 

= -2Llog2 - 2M 2 [log 2 - ±((3) A 2 - J + J'} +2J' + J- A 3 C(3)^ + 0(~). (B9) 

L * J 7T L 

It thus follows that the Fermi momentum A is defined by 

J' = J-log2 + IC(3)A 2 , (BIO) 
and that the A^ 4 ) become massless on this line at the point 

2J' + J = A 3 C(3)^. (Bll) 







We thus conclude that this phase boundary ends in the point 

(J, J') = (| log 2, -i log 2) 

parallel to the line 2 J' + J = 0. 



(B12) 
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